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Abstract-A novel design method for photovoltaic applications 
(PVA) is proposed based on Extreme Value Theory under climate 
change. Traditionally, one day or one year's climate data are 
used in the design of PV applications. However, the fact is that 
weather conditions in two years won't be the same.Therefore, the 
PV applications designed by the traditional way may fail in the 
extreme weather conditions. This paper illustrates a novel 
method to design a PV application reliable under the climate- 
change situation. 



Keywords-Climate Change; Extreme 
Application; Sizing; Climate Cycle 



Value Theory; PV 



1. The system obeys energy conservation law. 

2. The load is constant. 

3. The system is optimized by the vendor and the average 
power of the PV panel is used in calculation. 

4. The batteries are initially full. 

5. The batteries are only charged by the PV system at day. 

6. The load is supplied partly by the PV system at day and 
totally by the batteries at night. 



I. INTRODUCTION 

Design of Photovoltaic applications (PVA) is an active 
research area in the literature. A very comprehensive review 
about this topic can be referred to [1], The technical 
introduction of PVA can be found in [2]. The traditional 
methods of power system reliability evaluation are discussed in 
[3], Gordon [4] gave an optimal sizing technique based on the 
stochastic process. Egido and Lorenzo [5] presented an early 
review and proposed a method to take the advantages from 
both numerical and analytic approaches. Markvart et al. [6] 
proposed a sizing technique based on the extreme events. 

The design of PV applications mainly involves determining 
the number of PV plates and the number of batteries in the 
application. Most of them only use one day or one year's data 
to do the design task. However, the fact is that the weather 
conditions in a single site in two years will not be the same. 
Therefore, the PV applications designed by the traditional way 
may fail in the extreme weather conditions. 

Extreme Value Theory (EVT) [7, 8] is used in the modeling 
of extreme events occurred in a time series. Since the maximal 
climate cycle in a year normally cannot be predicted, it is taken 
as a random event and can be modeled by EVT. EVT has 
various findings that fit different situations. This paper mainly 
illustrates the design of PV applications in terms of block 
maximum. The realistic data from a real site are employed to 
the illustration. A comprehensive version of this approach can 
be referred to [9].The remainder of the work is described as 
follows: The assumptions for illustrating the approach are 
presented in Section 2. The mathematic preliminaries are 
addressed in Section 3. The illustration of the proposed method 
with realistic data is in Section 4. Finally, Section 5 is the 
conclusion. 

II. ASSUMPTIONS 

To simplify the illustration, it is assumed to satisfy the 
following assumptions: 



III. PRELIMINARIES 

A. Extreme Value Theory 

EVT is a method used for the evaluation of LPP X under 
observed data. The extreme value theory gives the following 
important results [10]: 

Theorem 3.1. Suppose X lt X 2 , ..., X n are iid random 
variables. If there are constants a n e R, b n > 0, M„ (the order 
statistics for X„) and some non-degenerate limit distribution H 
such that 

M -a 

limPr( — < x) = H(x) 

b n 

then H is one of the following extreme value distributions: 

f 0, x < 



Frechet: <t> (x) 



Weibull: v P (x) 



,« > 0, 



^e x > 



x < 



,a > 0, 



x > 



Gumbel : A a (x) = e , x e R 

A more general form called the generalized extreme value 
distribution (GEVD) of the three distributions can be expressed 



H 



x - /; 



_ (1+ ^ r i/f 



# "0, 



1 = 0, 



where £, is the tail index, ju is the position parameter, a is the 
scaling parameter. 

Theorem 3.1 states that the tail of any iid random variable 
with unknown distribution will asymptotic converge to the 
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GEVD. This is a very important result of precisely estimating 
LPP X mathematically. 

B. Calculation ofLPP x 
LPP X is defined as: 



LPP =Pr(M >x)=l-H c (—^-). 



C Climate Cycle Creation 

The climate cycle can be defined as the specific charge time 
to charge full the batteries. The following flow chart is 
proposed to create such climate cycle: 
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Fig. 1 The flowchart of creating climate cycles 

We define a symbol £ to do the summation. By definition, 
a climate cycle c, may start from the s, percentage of the PSH at 
the last day of the previous cycle c^ and end at the (1 - s i+ i) 
percentage of the PSH at the last day of this cycle. Let Tc t be 
the cycle time (i.e., the discharge time) of cycle c„ then it is 
expressed as 



I 24 = Tc 



s.e rn +12- 



"C-Dr 



- + 24(R-1) 



+(1-S i+ X+12- 



24(R-l) + (l-s, +1 )e . +12- 



1<!<Q, 



z=l. 



In unit of a day, it is denoted as 



Td, 



e ( ,_ 1)R , (S f -1/ 2) / 24 + R + e lR (1 / 2 - s i+1 ) / 24, 1< i < Q, 



R-l/2 + e t (l/2-s i+1 )/24, 



i = l. 



Tdf will be used in the later computation of LPP m 



Suppose there are Y years. According to the extreme value 
theory for block maxima, we define the random variables as 



Xj = max{Td,|c, is in year j}, 1 < j < Y. 



D. 



Sizing Formulas 

Let Tc max - max{ Tc, | 1 < i < Q}. How to choose a suitable 
charging time T and get the corresponding Tc max depends on the 
given weather condition. Empirical study shows that the 
relationship between Tc max /T and T is not linear, actually it is a 
nearly convex function. Figure 2 shows this relationship. 



g 
o 




Charging time (hrs) 
Fig. 2 The convex function realtionship 

Two sizing formulas are summarized from the observation. 

For PV panels: 



Mi- 



1 Tc 

-(— 
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For batteries: 



VI 



Tc 



RTE 



-(" 



-1). 



where RTE is round trip efficiency, P is the power of PV 
panels, VI is the power of batteries. 

IV. EXAMPLES 

Two sets of the solar radiation data in Florida, USA are 
illustrated: one is for the in-sample test measuring daily solar 
radiation for 30 years from 1961 to 1990; the other is for the 
out-of -sample test measuring daily solar radiation for 15 years 
from 1991 to 2005. Both data are obtained from the Renewable 
Resource Data Center of the National Renewable Energy 
Laboratory, USA. Assume that there is a constant daily load of 
4 kWh (48 V) uniformly distributed in a day to be powered by 
the PV system (205 W, 48 V for one plate) with batteries 
regulated, where 205W is the average power of the PV panels 
during PSH. Then, each module of the battery bank needs 4 
batteries (lead-acid type, 12 V, lOOAh) cascaded. The RTE of 
the battery is 0.9. 
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At first, the empirical analysis of the ratio Tc max l T vs. T is 
conducted. The minimal ratio 9.42851 is at 58.7 h (=T*) of the 
turn in this analysis. The corresponding length of maximal 
climate cycle is 23.06058 days Td * max . 

In in-sample test, we got the optimal PV panels are 9 plates: 

000 / 
24 [' 



Po = t (i+^(^-i))~(i+~- y 



= 1727.5Q2(W] = 8.43 plates. 
The optimal batteries are 4 modules: 
4000 



W = P D - i-o = 1727.502 



24 



1560.835(W) = 3.25 modules. 



The reliability estimation is 0.02981 which comes from the 
EVT fitting by MLE method. The fitting parameters are: \ - 
0.000001, ju = 1.048 and o= 20.947. Figure 3. is the fitted 




Fig. 3 The fitted extreme climate model 

This result is as good as that by the traditional sizing curve 
method, but with more mathematic insights. 

In out-of-sample test, we verify what happened to this 
design of PV applications. The climate statistics for the next 15 
years are in Figure 4. 
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Fig. 4 The extreme events in the next 15 years 

The results show that the optimal design is feasible for the 
future climates. However, three extremes are discovered by this 
method, which would occur at a frequency of per 317.5611, 



48.146 and 47.7783 years respectively. These are the evidences 
of climate change. 

V. CONCLUSIONS 

This paper illustrates the way of how to use the modern 
mathematics in the real applications. A novel design method 
for PV applications is proposed. Many interesting findings are 
obtained, which are never explored before. More realistic 
climate data are encouraged for the future verification of this 
method. 
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